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A Duffing oscillator in a certain parameter range shows period-doubling that shares the same 
Feigenbaum ratio with the logistic map, which is an important issue in the universality in 
chaos. In this paper a globally coupled lattice of Duffing flows (GCFL), which is a natural 
extension of the globally coupled logistic map lattice (GCML), is constructed. It is observed 
that GCFL inherits various intriguing properties of GCML and that universality at the level 
of elements is thus lifted to that of systems. Phase diagrams of GCFL are determined, which 
are essentially the same with those of GCML. Similar to the two-clustered periodic attractor of 
GCML, the GCFL two-clustered attractor exhibits a successive period-doubling with an increase 
of population imbalance between the clusters. A non-trivial distinction between the GCML and 
GCFL attractors that originates from the symmetry in the Duffing equation is investigated in detail. 



I. INTRODUCTION 

A globally coupled map lattice (GCML) is a system 
of generally a large number (N) of maps that iteratively 
evolve in discrete time (n) and interact via their mean 
field. One typical example is a homogeneous GCML of 
logistic maps [H-Q: 

N 

x.{n + !) = (!- e)fa{x.,{n)) + ^ ^ fa{xj{n)) 

This evolution can be regarded as an iteration of a two- 
step process. The first step involves parallel mapping 

xr'' = fix.,{n)),t = l,---,N (2) 

where randomness is introduced in the system if the non- 
hnearity parameter (a) of the map is set to be large. The 
second step involves interaction with a coupling constant 
e via the mean field: 

x^{n + 1) = (1 - eX"'* + eh{n), i = l,--- ,N 

Mn)^-^^™'^. (3) 

1=1 

Here the distribution of maps is contracted to the mean 
field h{n) at the rate of 1 — e by similarity transfor- 
mation, and all the maps arc forced to respect their 
average. Hence, by this averaging step, coherence is 
introduced into the system. Note that the relation 
X^iLi Xi{n + 1)/N — h{n) holds and the average is kept 
invariant at the interaction step in this model. This is 
important for the stability of the model. This GCML is 
simple but it shows us how the confiict between random- 
ness and coherence produces various interesting phases. 
(Kaneko [l|). 

However, in a way, GCML is oversimplified. The real 
physical process (say, pattern recognition in the brain) 



occurs in continuous time. At this point, it is tempting 
to construct a concise model of the globally coupled fiows 
(GCFL) and study to what extent it inherits various in- 
triguing properties of GCML. This study aims to search 
universality in complex systems ( GCML and GCFL ). 

How to choose the flow (as an element of GCFL) that 
replaces the logistic map (as an element of GCML)? Ac- 
cording to the Poincare-Bendixon theorem, we must con- 
sider a three- or higher-dimensional flow as an element 
to exhibit a chaotic behavior at the element level if the 
flow is autonomous, or alternatively, one may consider a 
two-dimensional flow under an external force. In this pa- 
per, we choose, for simplicity, a two-dimensional Duffing 
oscillator under an external periodic force 

dx dy , Q , , , ^ 

-=,, - = -ky-x^-,Acost (4) 

with A = 7.50 Indeed, it is known that the Duff- 
ing flow shares the same period-doubling bifurcation to 
chaos with logistic maps, sharing the same Feigenbaum 
ratio 4.6692016- • • 0,0]. The reason for this univer- 
sality is as follows: the two-dimensional flow of this 
model produces a one-dimensional iterated map on the 
Poincare section and, therefore, the bifurcation of this 
model is in one-to-one correspondence with the bifurca- 
tion of the one-dimensional iterated map. Besides, all the 
one-dimensional maps (with a singly peaked function) are 
governed by a universal function at the limiting bifurca- 
tion. The above scenario is crucial in the universality in 
chaos 0. In this paper, we numerically study to what 



This simple parametrization is studied by Yoshisukc Ueda [J, 
who found chaos in nonlinear oscillators by analog computer sim- 
ulation. In addition, Feigenbaum discusses the Duffing oscillator 
in this form . We have verified that our findings are essentially 
independent of the choice of the potential shape. The double- 
well potential Duffing model also follows the symmetry JSj, and 
similar to our single-well Duffing GCFL, the double- well Duffing 
GCFL also exhibits two-clustered attractor under the symmetry 

m ■ 
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extent the universality at the element level (the univer- 
sality between the logistic map and the DufRng oscillator 
with respect to the bifurcation to chaos) extends to the 
similarity between the systems, i.e., the logistic GCML 
and the Duffing GCFL. 

In fact, the above universality (at the element level) ex- 
tends to a certain higher-dimensional flow and relates it 
down to a one-dimensional map. For instance, the three- 
dimensional Lorenz model Q (with ct = 10 and 6 = 8/3) 
has three prominent r bands, in which it shows period 
doubling Q . The key to this issue is that the Lorenz sys- 
tem is dissipative. Hence, the effective dimension of the 
orbit reduces to two in this regime, and on the Poincare 
section, it is equivalent to a one-dimensional iterated map 
0, 0]. Our preliminary analysis of the coupled Lorenz 
system shows intriguing results: (1) For strong coupling, 
a two-clustered attractor is formed just as in the Duffing 
GCFL. (2) For weaker coupling, the distribution of flows 
(at any time after the transient process) forms a one- 
dimensional closed string (single type) in the three- 
dimensional phase space, and as we reduce the coupling 
further, the string bifurcates at the critical point into 
two closed strings (two strings linking each other). 
This type of topology change is known to occur for a sin- 
gle Lorenz attractor [l^l ■ We defer the discussion of the 
Lorenz GCFL to a forthcoming paper. 

Note that the Duffing equation ^ has symmetry un- 
der the transformation 



differs from most of the other synchronization models of 
flows, where elements interact via the time derivatives 
{xi (t) ,yi{t),- ■ ■) • The manner in which the non- 

linearity of both GCFL and GCML is matched is de- 
scribed in detail. Then, we study the phase structure of 
the Duffing GCFL by varying the coupling e as a free pa- 
rameter. We find that the phase diagram of the Duffing 
GCFL resembles that of the logistic GCML. The univer- 
sality at the level of elements extends to that of coupled 
systems. In Sec. lIIIl we focus on the two-clustered regime 
of the Duffing GCFL. Here flows divide themselves into 
two synchronizing clusters and the orbits of the two clus- 
ters bifurcate as the population imbalance between the 
two clusters increases. This bifurcation precisely corre- 
sponds to the phenomenon in the two-clustered regime 
of the logistic GCML discovered by This may also be 
considered as an extended universality. However, there is 
a subtlety owing to the attractor pairing. We will clarify 
this issue in Sec. lIIIR . (Fig. Only after considering 

this distinction can we fully understand both similarity 
and dissimilarity of nonlincarity reduction in GCML and 
GCFL, as discussed in Sec. 11113 (Fig. [H] and [TOl) . 



II. DUFFING GCFL 



A. Construction of Duffing GCFL 
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, which allows for pairwise attractors. 
{xi{t),yi{t)) is a solution of (|4]), then 



(5) 

That is, if 



(X2(i),y2(i)) = (r,®P)(xi(t),yi(i)) 

= (-xi(< + ^),-yi(t + ^)), (6) 

is also a solution, where the symbols Tj^ and P, respec- 
tively, denote the time translation for half period of the 
external force and the parity transformation in the two- 
dimensional phase space. Depending on the initial value, 
a Duffing flow may fall in one of the two attractors. On 
the other hand, the logistic map has only a single attrac- 
tor independent of the initial point. This difference at 
the element level is reflected in the final structure of the 
attractors of coupled models and we will discuss this in 
detail. 

Some results in this paper were reported at AROB 
13 [ll|. This paper presents our completed work with 
a deepened understanding of GCML-GCFL correspon- 
dence (especially relation pT|) and (HU). On the basis 
of high statistics simulations, we report below detailed 
phase diagrams of the Duffing GCFL and the bifurcation 
diagram with respect to the cluster composition ratio ??. 

The rest of this paper is organized as follows. In Secllll 
we first construct a Duffing GCFL along the same lines 
as GCML. We put Duffing flows in all-to-all interaction 
via their mean field. This interaction acts on the "coor- 
dinates" of fiows {xi{t) , yi{t)) . At this point, this model 



We construct the Duffing GCFL by replacing the maps 
in GCML by Duffing flows. The flrst step © now be- 
comes 



xr''(t + At)=x,{t)+y,{t)At 

yr^it + At) = y,{t) + [~ky,{t) ~ {x,{t)f + Acos{t) 
z = l,--- ,iV 

and the second ([3]) becomes 

x,{t + At) = (1 - eD)xr''it + At) + eoh,. 

y,{t + At) = (1 - eD)yr\t + At) + eohy, i = 1, • • 



At 
(7) 



1 ^ 



At), 



hy = 



1 " 



y™'*(i + At). 



(8) 

The interaction in this step is via the mean field of the 
system and the change in the coordinates is directly cal- 
culated by a similarity contraction, similar to that in 
GCML. Hence, let us designate this model as GCFL. 
A similar model was investigated in the earliest stage of 
the chaos synchronization study [iJilili (Eq.(5) in jisj). 
There, the synchronization of the chaotic elements was 
the main concern, while in this study, our interest lies in 
how various synchronizing phases appear from the origi- 
nally chaotic fiows via nonlincarity reduction due to av- 
eraging interaction. Interesting studies have also been 
conducted on phase synchronization (the amplitude may 
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be left unsynchronized) [l^ [31 . Here we study clus- 
ter formation (synchronization in both phase and ampli- 
tude). An extensive account on the synchronization can 
be viewed in [l3j . 

To investigate the phase structure of the Dufhng 
GCFL, we need some rough guide of the regions to ex- 
plore the vast (k, e) parameter space. 

Let us first discuss the case of nonlinearity parameters. 
To study nonlinearity reduction (i.e., the battle between 
chaos and coherence), the nonlinearity parameter of the 
elements should be set high in the random chaotic regime. 
For a logistic map, the threshold to chaos is ath = 1-401, 
so let us choose a = 1.7, 1.8, 1.9 as typical values. As for 
a Duffing flow, note that the parameter fc, the friction 
constant, is an anfi-nonlinear parameter. From univer- 
sality in chaos in the sense that two elements in the 
same universality class have the same bifurcation trees, 
a natural choice for comparison is to align the bifurca- 
tion tree of a Duffing flow with that of a logistic map, as 
shown in Fig. [TJ 

In fact, there is a subtlety in the Duffing flow. At very 
low nonlinearity (very large friction constant fc), the pe- 
riodic orbit is self-symmetric in the xy plane, and with 
an increase in nonlinearity, it changes the topology, split- 
ting into two periodic orbits {A+ and A_) at {k ~ 0.66), 
which are mirror symmetric of each other (see insets in 
Fig. [1]). We consider that each of the A± corresponds 
to a flxed point of a logistic map and align the Duffing 
tree with the logistic tree such that the second and third 
bifurcation points match between the two trees. 

It is interesting to test how the Lyapunov exponents 
compare between two models after matching the bifur- 
cation trees. The Duffing flow induces a Floquet map at 
the Poincare shot, which is taken at every period of the 
external force. Therefore, we expect a correspondence 



,x27r 



(9) 



The agreement between the exponents is remarkably 
good except for the details of windows, as observed in 

Fig.m 

Up to now, we have been considering matching at the 
level of the basic elements of two models. Now, let us 
consider the coupling ed- In GCML, the similarity con- 
traction by 1 — Emap occurs at each interaction step. On 
the other hand, in the Dufling GCFL (I?]),®, the con- 
traction by 1 — Ed occurs at each At, and At, in turn, 
has to be chosen sufiiciently small to guarantee flnite dif- 
ference approximation. (We typically use At = 10^^.) 
Therefore, the coupling ed must be chosen to be small 
enough such that the accumulation of the contraction 



^ We are aware that the universality holds at the bifurcation limits, 
so in principle, it would be better to consider the matching at the 
bifurcation points as high as possible. However, such matching 
using narrow strips is bound to produce large global errors. For 
a few different alignment choices, we have verified that the main 
body of our results is unchanged . 



effect in GCFL roughly amounts to the contraction in 
GCML. That is, very crudely 



1 -emap ~ (1 -Ed)^ , 



(10) 



where tc is a certain time scale of order one. If one may 
simply carry over the correspondence at the level of el- 
ements to that of the systems, one iteration step of a 
map corresponds to a 27r evolution of the flow so that 




k=0.39 



k=0.33 



k=0.29 



FIG. 1. (a) Bifurcation tree of a logistic map. (b) Bifurca- 
tion tree of a Duffing flow. For k > ko = 0.660, the Duff- 
ing attractor is unique and left-right symmetric (not shown), 
while for ko > k > ki — 0.380, it becomes bivalent (either 
A+ or A_), eacli being chiral (see insets). With a further 
decrease in k, each of A+ or A_ bifurcates. The threshold 



for A 



' — !> ' is denoted as k„. Arrowed dashed lines 
connecting (a) and (b) indicate that the second and third 
bifurcation points are chosen for the matching of a and k. 
Solid lines show A, B, and C points, which are used for the 
comparison of GCML and GCFL. 
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FIG. 2. Comparison of the element nonlinearity in terms of 
Lyapunov exponents. Upper: the exponent of a logistic map 
as a function of a in IT}. Lower: the maximum exponent of a 
Duffing flow as a function of fc in Q. Dashed lines indicate 
the second and third bifurcation points, which are chosen for 
the matching of a and k, and solid lines show A, B, and C 
points, which are used for the system comparison (cf. Fig.[TJ. 



one might assume tc = 27r. However, the systems arc 
evolving while interacting in nontrivial ways (the crucial 
point is the existence of a pairing attractor structure in 
the Duffing flow) and the above naive expectation is un- 
justified. 

We proceed below without any assumption on the 
value of the e (and tc)- We regard e as a free param- 
eter of the model and only use as a guide to explore 
the GCFL phase diagram. We will find below that in 
the Duffing GCFL, phases quite similar to those of the 
logistic GCML are realized and the proper choice of tc is 
as given in ([H]). 

Now, we are ready to investigate the Duffing GCFL. 
Fig. [3] shows a typical run in the Duffing GCFL with 
N = 100. Here the model (fc = 0.241 (case A) and e ^ 
1.1 X 10~^) is started at a random initial configuration, 
and after about 10 cycles of external periodic force, the 
flows divide themselves into two synchronizing periodic 
clusters. This is a typical run for the two-clustered regime 
of the Duffing GCFL. Let us now compare the overall 
phase structures between the two models. 



B. Phase diagrams 

1. GCML phase diagram 

The phase structure of GCML ([T|) is explored by 
Kaneko; refer to [l| for a detailed phase diagram in the 
(a, e) plane. To facilitate a comparison with the GCFL 
phases, we reproduce in Fig. U the GCML phases at A, 




FIG. 3. A typical run in the Dufling GCFL. N = 100, 
k = 0.241 (case A) and e = 1.1 x 10"'* {At = 10"^). Up- 
per: periodic external force F{t). Lower: synchronization of 
Dufling flows into two periodic clusters. Only y coordinates 
are shown. 




FIG. 4. GCML phase diagrams at points A, B, and C. 
(Kaneko With decreasing e, one observes three out- 

standing regimes (coherent chaos, periodic two-clustered, and 
turbulent regimes, respectively); in between these, there are 
periodic-clustered and partially ordered regimes. With in- 
creasing nonlinearity parameter a, a larger coupling e is re- 
quired to maintain the same dynamics. 



B, and C points. There are three outstanding phases: 

(1) The coherent chaos state at large e. The maps are 
strongly bunched together in one cluster {xi{rL) = 
X{n), i = I,-- - ,iV). Then, h{n) = X{n) and, 
therefore, the interaction ([3]) becomes immaterial. 
All the maps evolve in a bunch just as a single 
logistic map at the original a. 

(2) The two-clustered phase at intermediate e. The 
final maps divide into two clusters(^ = 1,2), 
(a;f (n) = X^^{n), i = 1, • • • ,iVp), and the two clus- 
ters, with populations A^i and N2, respectively, os- 
cillate periodically opposite in phase to each other, 
thus keeping the fluctuation of the mean field h{n) 



We have cross-checked these phase structures. 
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FIG. 5. Phase diagrams of the Duffing GCFL at A (A; = 0.241), B (A: = 0.228), and C (A; = 0.215), corresponding to GCML at 
A (a = 1.7) B (a = 1.8), and C (a = 1.9), respectively. (Cf. Figures [l] and ID) The time step for the finite time approximation 
is At = 10~^ and the coupling constant en is varied from to 24 x 10^^ with increments of 0.5 x 10~^. At each e, 50 random 
configurations are generated and the final attractors after t = 500(27r) are examined for cluster structures and orbit types. The 
insets exhibit the orbits of the flows in the final attractor in the respective phase. Basin volume ratios of the final attractors 
are exhibited by the partition of the bar width. The attractor of the type coherent- chaos at the strong coupling regime changes 
its topology. At the low nonlinearity side of the element flow (case A), the attractor is chiral, while at the high nonlinearity 
side (B and C), it is self-dual. With increasing nonlinearity (decreasing k), a larger coupling e is required to maintain the same 
dynamics. 



small. This opposite phase motion is a solution for 
stability. A remarkable property found by Kaneko 
[H is that the population ratio -d = A^i/(A^i + iV2) 
acts as a new bifurcation parameter. 

(3) The turbulent phase at very small e. In general, the 
number of clusters is proportional to N and maps 
evolve almost independently of the randomness of 
the original a. 

In addition, in a region marked as periodic clusters^ sev- 
eral clusters ((2,3), (2,3,4), (3,4),---) are formed with 



* However, statistical analysis reveals hidden coherence — a long- 
time-scale correlation 0] . Furthermore, there occur drastic peri- 
odic motions of maps in stable clusters at specific e values tuned 
with a's — periodicity manifestations such as p3c3 and p3c2 ll7l| . 
The reflection of these manifest cluster formations dominates all 
over the turbulent regime of GCML. These are due to the folia- 
tion of periodic window dynamics of the element maps l lIV A It . 
We refer to [H, where the foliation at the small e region is caught 
as a tongue-like structure. 



decreasing e, and in partially ordered phases, the basin 
volume is partitioned by the dynamics of two surrounding 
regions. 



2. GCFL phase diagram 

Fig. [5] shows the phase diagrams of the Duffing GCFL 
for three choices of the friction coefficient, A:=0.241, 0.228, 
and 0.215 corresponding to A, B, and C, respectively, 
in Fig. [T] In the figure, basin volume ratios of the 
final attractors are exhibited by the partition of the 
bar width. For instance, at fc = 0.241 (point A) and 
e = (10.4 - 13.4) X 10~^, the flows almost always (with 
10% exception) fall into two-clustered attractors (both 
periodic (P(8)P) as shown in the inset), while at the same 
k but at e = 15.5, the basin volume of the P(8)P attractor 
and that of the coherent chaotic attractor are almost the 
same (the center of a partially ordered phase). Compar- 
ing Fig. [S] with Fig. 01 we observe that the phase struc- 
ture of the Duffing GCFL is generally the same as that of 
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GCML. The three basic phases of GCML, namely (1) the 
coherent chaos, (2) the two-clustered phase, and (3) the 
turbulent phase, are also realized in the Duffing GCFL 
and in the same order as coupling strength. Moreover, 
the tendency that the regimes shift to the higher coupling 
with increasing nonlinearity is also the same. (Note that 
the lower friction k implies higher nonlinearity.) The 
agreement in the basic phase structure implies that the 
Duffing GCFL inherits the intriguing properties of the 
logistic GCML and that the universality at the element 
level may be lifted to that at the level of globally coupled 
systems. 

Below, we list the nontrivial properties of the Duffing 
GCFL that are not present in GCML: 

(1) Symmetry of the Duffing equation under ([5]) leads 
to an attractor consisting of pairwise clusters. The 
clusters may be called mirror pairs (but with a shift 
of TT in time.) (On the other hand, in the GCML 
two-clustered regime, clusters evolve in essentially 
the same orbit with a shift of one step.) 

(2) As an intriguing consequence of (1), we observe not 
only the P (g) P state but also the C (g) C and C P 
states. 

(3) If the coupling is high enough, the flows are 
bunched in a cluster and the cluster evolves with 
the original high nonlinearity of the elements (co- 
herent chaos). At this point, the case is the same 
as that for GCML. However, the Duffing flow at 
high nonlinearity shows two types of attractors. If 
the nonlinearity is not so high, the chaotic flow is 
confined in either the left or right polarized orbit 
(chiral type), while at the higher nonlinearity, the 
chaotic flow extends over the joint of the above two 
orbits (self-dual type). This is reflected in the co- 
herent chaos state. (Compare the coherent chaos 
in A with that in B and C). 

Hereafter, we focus on the interesting two-clustered 
phase, which is realized in both GCFL and Duffing 
GCFL. 



III. TWO-CLUSTERED DYNAMICS OF THE 
DUFFING GCFL 

A. Attractor t?-bifurcation 

1. -bifurcation in GCML 

In the two-clustered regime, the GCML attractor con- 
sists of two clusters mutually oscillating opposite in 
phase. Let iV+ be the number of maps in the positive clus- 
ter (the one consisting of maps with Xi{n) > h(n) at even 
iteration step n) . Kaneko found that with an increase (or 
decrease) in the population ratio = N^/{Njf.-\-N-), the 
cluster attractor repeats successive bifurcations until the 




6=0.23 9=0.5 6=0.77 

type B type A 



FIG. 6. Two-clustered attractors in GCFL with iV = 100, 
k — 0.241, e = 11 X 10"^. The cluster population is examined 
att = 100 x 27r) to determine 6* = N+/{N++N-). The motion 
of majority (minority) cluster thereafter is depicted by white 
(black) circles on its orbits depicted by a bold (dashed) line 
at every 7r/2 step. In type A (B) attractor with & > 0.5 
(i9 < 0.5), the positive (negative) cluster is the majority. One 
finds an exact symmetry pip and an approximate symmetry 



imbalance reaches a certain threshold beyond which the 
system goes into chaotic transients [11, 3 . Let us call this 
phenomenon iJ-bifurcation. It implies that the cluster at- 
tractor is labeled by three parameters (a, e, and and 
it may lead to a possible new approach of information 
processing. 

2. Attractor symmetry and ^-bifurcation in GCFL 

Let us investigate if our Duffing GCFL inherits this 
intriguing property. The answer is yes, but, in order to 
fully understand the i9-bifurcation in GCFL, we have to 
take account of an attractor symmetry coming from the 
symmetry ([5]) of the Duffing equation. 

In Fig. \6\ we show three typical events in the GCFL 
two-clustered regime. The central event is basic. It is 
1? = 0.5 and is of period one. If a time translation of tt is 
applied on one of the two clusters, its position becomes 
exactly mirror symmetric to that of the other cluster. 
The other two events are d = 0.77 and i? = 0.23. The 
cluster orbits are indeed bifurcated and are of period two 
[T = 47r). Let us call the attractors type A and type B 
if ?? > 0.5 and d < 0.5, respectively, and two attractors 
as partner attractors ii i!)a + = 1- The attractors 
shown at the right and the left of the basic attractor, 
respectively, are thus partner attractors. Inspecting the 
motion of clusters in Fig. [6l we flnd an interesting sym- 
metry that holds between the partner attractors: 

(T, ® P)(XA(t),YA(t))[^l = i~XAit + tt), -YA{t + ^))(<^) 
= {XB{t),YB{t)Y^'^ {C = L,S), (11) 

where C ~ L, S denote the large (majority) and small 
(minority) clusters, respectively. This symmetry holds 
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for majority and minority clusters separately. Note that 
at the limit — > 0.5, this relation reduces to the symme- 
try observed in the central event. This, in turn, implies 
that we also have an approximate symmetry relation 

« {XB{t),YB{t)Y^^''\ (12) 

as confirmed by Fig. |6l 

Given that the dynamics of the N flows is reduced by 
synchronization to that of two clusters, that is, assuming 
the two-clustered configuration, it can be seen that the 
symmetry ([6]) of the single Duffing oscillator is lifted to 
(fTTj) . We show below that if ([TT|) holds at time t, then it 
also holds at time t + At. 

The GCFL evolution is an iteration of the two-step 
process, so let us start with the first step. For the two- 
clustered configuration, this leads to 

X'^id^^ At) = F{X{t),Y{t)) = X{t) + Y{t)At, 
Y"'"^{t + At)=GiX{t),Y{t),t) (13) 
= Y{t) + [-kY{t) - X{tf + A cost] At. 

This applies to both majority and minority clusters and 
to both type A and B attractors. First, let us examine 
yrmd^ where G has an explicit t dependence. We have 

Yr\t + ^t) = G{XA{t).YA{t),t) 

^G{-XB{t + n),-YB{t + ^),t) 

= ~G{XB{t + ^),YB{t + TT),t + TT) 

= -y™'^(i + 7r + At). 

The first and forth equalities are from (fT3|) . the second 
follows from the assumption at t, and the third from the 
symmetry of the Duffing equation ([6]) used in the first 
step. Thus, we obtain lX"^(t-h At) = -Y^^'^[t + TT+ At). 
For X""'', F does not involve t and we immediately find 
that XJ^-^it + At) = -X^''^{t + n + At). Thus, we have 
shown that m]) holds for both X™'^ and F™*'*. The 
next step is the interaction through the mean field. For 
the two-clustered configuration, for the X coordinates, 
we have 

X^^\t + At) = (1 - £_D7?s)^^-^'™''^ + eDt?5^^^'™''^ 

X(^'(t + At) = (1 - sd^l)X^'^-"'"^^ + ed^lX''^''^"^^ 

and the same also holds for the Y coordinates. As these 
are a linear mapping, (|lip is apparently kept. Thus, we 
have shown that (fTTI) is the symmetry of the GCFL evo- 
lution. 



3. A comparison between the GCML and GCFL 
^-bifurcations 

With this symmetry in mind, let us compare the "d- 
bifurcation in GCML and that in GCFL in detail in 
Fig. [71 We use a representation devised by Kaneko [l| . 
As for GCML, we draw all values of the maps at d in the 
final cluster attractor at all even time steps. First, con- 
sider an event of a two-clustered attractor with d — Q.h 



(the evolution is not drawn). This is the basic attrac- 
tor in the GCML two-cluster regime, in the sense that 
it is free from i?-bifurcation. The two clusters oscillate 
mutually opposite in phase around their mean field, and 
each one returns to its previous position at every two 
steps. Thus, the basic motion produces two points at 
i9 = 0.5. In Fig. El (a) a sample event with ?9 = 0.41 
is exhibited. The orbits arc bifurcated to period four 
for high population imbalance, and they contribute four 
points (two for each cluster) to the iJ-bifurcation diagram 
[(b)]. The omission of odd step data serves to display the 
i9-bifurcation of two clusters separately in one diagram. 

As for the Duffing GCFL, the basic motion at d — 0.5 
is period one (T = 27r), as discussed in Fig.|6l Therefore, 
we take the Poincare shot of all flows at every multi- 
ple of 27r and plot them at the d to construct the d- 
bifurcation diagram. In Fig. |71^c) and (e), we show a 
sample event with = 0.23. Due to the population im- 
balance, the two-clustered attractor is bifurcated to pe- 
riod two [T — Att) , and the cluster orbits again contribute 
four points in the ??-bifurcation diagram. In GCFL, the 
mean field hx^y{t) oscillates because of the external peri- 
odic force and the two clusters mutually oscillate around 
this oscillating hx^y{t). However, at the Poincare shot, we 
clearly observe that the two-clustered attractor in GCFL 
shows a remarkable 79-bifurcation. which is quite similar 
to that in GCML. 

It is worth summarizing two clear distinctions between 
the GCML and GCFL ?9-bifurcation. 

(i) The basic motion at i? = 0.5 in the two-clustered 
attractor regime is period two in GCML and period one 
(T = 2n) in GCFL. That is, two iteration steps in GCML 
correspond to one rotation in the Duffing GCFL: 

2-clustered regime 
GCML ^ Dufiing GCFL (14) 

n n + 2 t t + 27r. 

The "iJ-bifurcation occurs successively from this basic at- 
tractor with increasing population imbalance. 

The matching of the nonlinearity between the elements 
is fixed in III Al on the basis that one iteration of a free 
map corresponds to one rotation (T = 2it) of a free Duff- 
ing fiow. There is no contradiction because in (|T4l) . the 
correspondence between two systems is considered, where 
interaction is in action and system nonlinearities are re- 
duced. We discuss this point in detail in the following 
subsection. 

(ii) In GCML, each of the two clusters evolves in its one- 
dimensional orbit and two orbits of the clusters approxi- 
mately agree with each other (precisely if ?9 = 0.5) mod- 
ulo a shift of one step. In the Duffing GCFL, on the other 
hand, the cluster attractor is two dimensional under the 
T (S) Ptv symmetry pT|) . Both (i) and (ii) together real- 
ize quite similar ?9-bifurcation diagrams for GCML and 
GCFL. 
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FIG. 7. The i^-bifurcation diagrams (right) associated with the evolution plots of the sample events (left) illustrate how diagrams 
are calculated. The logistic GCML [a, b] and the DufHng GCFL x-coordinates [c,d], y-coordinates [e,f]. In the evolution plot, 
majority and minority cluster orbits are depicted by solid and dashed lines, respectively. Their positions, at every even step for 
GCML and at every even multiple of tt for GCFL, are depicted by white and black circles for the majority and minority clusters, 
respectively, and they contribute to the i9- bifurcation diagrams altogether four points (two for each cluster), as indicated by 
horizontal fine dotted lines. 



4- Polarization of the ■& -bifurcation tree in GCFL 

Let us add a brief remark here. The GCFL 
bifurcation diagram is strongly polarized (the larger d 
forks arc missing) in contrast to the unpolarized GCML 
diagram. This is possible since the external oscillatory 



force may favor type-B rather than type-A attractors in 
the formation process, while GCML is autonomous. In 
order to overcome unwelcome missing bifurcation trees 
in the larger -d side, we can use the relation ([TT]). If we 
measure the type-B event (with some = i?o) at every 
odd multiple of tt, it supplies equivalent data for the type- 
A event (with d — 1 — z?o) a-t every even multiple of tt. 
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FIG. 8. GCFL i9- bifurcation diagrams, (a) is the same with Fig. [7(d) (flows at every even multiple of 2ii) except that the 
number of initial configurations is limited to 10'' to show i9 polarization. In the event formation probability below, one observes 
that type B is clearly more frequently produced, (b) is deduced via ([5]) from the same 10'' data, but orbits at odd multiples 
of 2tt axe converted by the relation (|11|) . (c) is our final GCFL i9-bifurcation tree constructed from a total of 6.4 x 10^ events; 
both even and odd steps data are combined). 



Fig. [5] presents the GCFL polarization- free i9-bifurcation 5. A comment on a posi-nega switch 

diagram thus obtained. To our knowledge, this is the 
first observation of i?-bifurcation in the globally coupled 
Duffing flows and its comparison with i^-bifurcation in 

GCML '^^ would be interesting to investigate if a posi-nega 

switch found in the GCML [l| can also be realized in 
GCFL, but we defer this for the following reasons. In 
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GCML, the two periodic clusters have negative Lyapunov 

exponents!! In several steps of iteration, the variance 
of map coordinates in each cluster becomes less than 
the machine epsilon One could introduce low level 
noises to avoid numerical collisions, but then the analysis 
of the posi-nega switch will involve technical subtleties. 
The GCFL attractors have also negative Lyapunov ex- 
ponents and share the same difficulty in the numerical 
analysis. Furthermore, we have checked that the GCFL 
two-clustered attractors are robust in the sense that the 
threshold value "dth to the decay into a chaotic transient 
process is quite high {"ffth » 0.9) (or quite low dth ~ 0.1). 
This is a negative indication for an implementation of 
the switch, because, in order to enhance the population 
imbalance beyond the threshold, many elements must be 
transferred from the minority to the majority cluster and, 
then, the above collision problem becomes more severe in 
the GCFL than in the GCML. 



IV. TWO-CLUSTERED DYNAMICS OF THE 
DUFFING GCFL 

A. Nonlinearity reduction 

The periodic attractor in the two-clustered regime is 
a consequence of the reduction in the high nonlinear- 
ity of the elements caused by the averaging interaction 
via the mean field. We first recapitulate this issue in 
GCML, where some quantitative understanding is possi- 
ble on the basis of the Perez-Cerdeira transformation for 
the quadratic map. We then numerically investigate the 
case of GCFL in some detail. 



1. Case of GCML; foliation curves of periodic window 
dynamics of element maps into the (a, e) plane 

In GCML, all element maps have high nonlinearity, 
and the maps would evolve independently in random mo- 
tion if the interaction would be switched off. However, at 
every iteration step, the map coordinates are uniformly 
contracted to the mean field and coherence is introduced. 
The nonlinearity of the system is effectively reduced and 
the periodic clustered attractors are formed. 

Now, interestingly, there is a case in which this reduc- 
tion can be quantitatively estimated. This is the case 
of the maximally symmetric cluster attractors (MSCAs) 



formed in the turbulent regime. We refer to [l7[ for a 
detailed account. In the MSG A, the maps synchronize in 
almost equally populated c clusters that oscillate in pe- 
riod p = c around their mean field. The most prominent 
one is the p3c3 MSCA. In such an MSCA, the mean field 
is almost constant because of the population symmetry, 
and hence, it is possible to estimate to what extent the 
nonlinearity of the element maps is reduced. In fact, as 
the mean field is approximately a constant hn ~ h*, the 
GCML evolution equation leads to 

Xi{n+l) = (1 -e) [1 -a(x,(n))2] + sh* , i = l,--- ,N. 

(15) 

Then, the scaled maps defined by the Perez-Cerdeira 
transformation 1231 



Xiin) 



(16) 



\~e + eh* 
all equally follow the same scaled map 

y,{n + l) = l-b[x,{n))\ i^l,---,N. (17) 
The nonlinearity is reduced from a to b = ra, where 

r = (l-e)[l-e(l-/i*)] < 1 (18) 

is the rate of nonlinearity reduction. A good rule of 
thumb is r « 1 — 2e. 

We can go a step further. Taking the average of ([2 
over the elements, we obtain 



h* 



1 - e + eh' 



= (y) = vl 



(19) 



, where yl is the long-time average of a single map yb(n) 
with nonlinearity b. Solving ([T5|) and (|19p for a and e 
at a given b and r, one obtains the parametric form of a 
curve (labeled by b and develops with r) that describes 
the foliation of dynamics of a period p window into the 
(a, e) plane. Q 



a('')(r) = ^, 
r 



e^''\r) = 1 



^r{l-yl) + 



(20) 
(21) 



We generally observe the period p MSCA (r > 0.95) and 
its remnant state {r < 0.95) in the (a, e) region, which lie 
between the two foliation curves emanating from bi and 



^ Precisely, for the two-clustered attractor state with populations 
A^i and N2, the Lyapunov exponents generally divide into two 
sets. The first one consists of two exponents, A^i — 1 and N2 — 1 
fold degenerate, respectively. These are responsible for the con- 
traction of the maps to the center of the respective cluster. For 
the two clusters of concern, we have measured these as negative 
l20t . The second consists of two exponents that are responsible 
for the stability of the cluster orbits and they are also negative. 



^ (y) is, for the MSCA, an equal-weight average of the clus- 
ter orbits over p = c clusters. But, again for the MSCA, 
this is nothing but the average of the cluster orbit points 
Yj{n) over period p. This, in turn, is also the case with y^. 
Namely, (y) = Ei=i '^/{«)(Vri) = JLltV;-' Yi{n){VI) = yl = 

^ Om formulation of the foliation curves [17| was inspired by l2lll . 



11 



Logistic GCML 



Duffing GCFL 



x10-' 



0.5 



0.4 



Coherent 




20 



1 n n \ 



Z3 




fl] H 



1.4 1.5 1.6 1.7 1.8 1.9 2.0 



J L 



non-linearity 



0.20 0.21 0.22 0.23 0.24 0.25 0.26 0.27 

< k 

non-linearity 



FIG. 9. The band of the two-clustered dynamics of GCML (left) and GCFL (right) in the parameter spaces (a,eniap) and 
(kjEo), respectively. (For comparison, the e_D-axis is adjusted to the Smap-axis using the relation (|22|) .) In the black band, 
the formation probability of the two-clustered periodic attractor is more than 90while in the adjacent gray regions, the basin 
volume is shared more by nearby dynamics (coherent chaos in Zi and Z^, multiclustered attractors in Z2,C ® C in Z4,), and 
toward the outer edges of the gray bands, the partition of the two-clustered periodic attractor decreases down to 10% (Fig. [Sj ■ 
The judgment of final states is performed by inspection after 10'' iteration steps from an initial random start. If it is performed 
after 10^ steps, Zi is partially ordered. For GCML, the foliation curves of the dominant windows are also shown; in particular, 
61 and 62, respectively, emanate from the first and second bifurcation points (pi — )■ p2 and p2 — >■ p4) (ref. [itI. |20||'|. 



&2 (indicating nonlinearity at the p window's opening and 
at the first bifurcation in the p window, respectively). 
Ifl In the adjacent region (with larger e), period p, c- 
clustered states (with p > c) are formed. The foliation 
curves are exhibited in Fig. |9]-GCML. 

Wc should now clarify the logic. The foliation curves 
are derived assuming the existence of the MSCAs. If 
a p = c MSCA is formed, it must be within the (a, e) 



region dictated by the foliation curves of the period p 
window. Furthermore, we can prove that at the MSCA, 
all the Lyapunov exponents are negative and the MSCA 
is a stable attractor (T^]. However, the converse is not 
necessarily true. There is no guarantee that the MSCA is 
formed everywhere in the range specified by the foliation 
curves; indeed for r < 0.95, we observe only the remnants 
of MSCAs {p > c clustered attractors) as the valleys 
(peaks) in the distribution of the mean squared deviation 
of the mean field h{n) [l3,[ll- 



The most prominent p3c3 is an exception. It is observed down 
to r ~ 0.84. 



Now, in the two-clustered regime of concern, the maps 
divide themselves into two clusters, and the two clus- 
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ters oscillate mutually around the mean field. Here the 
attractor is robust in that it allows the population imbal- 
ance accompanied by the attractor bifurcation. However, 
in the particular case oi d = 0.5, we can regard it just as 
a p = c = 2 MSCA. With this concept, we can apply the 
above formulation to predict the possible region in which 
the two-clustcrcd attractor may be formed. 

In Table. U we compare the observed e range of the 
two-clustered regime with the predicted bounds from the 
foliation of the period-two window dynamics. For the 
latter, the two curves from bi = 0.75 and 62 = 1.25 
are used, which are the opening point of the p2 window 
(the threshold of pi — ?• p2) and the threshold of p2 
p4, respectively. Table. |I] shows that the observed range 
{£2, £3) is inside the predicted bounds. Moreover, Fig. [9]- 
GCML shows that the observed band for the formation 
of two-clustered attractors runs just in the center of the 
above two foliation curves. 



TABLE I. Two-clustered regime in GCML 



Nonlincarity'^ 


observed 


£-range^^ 


Prediet ion (bounds^ 




ei - £2 


- £3 - £4 




A (a=1.70) 


0.16 - 0.19 


- 0.25 - 0.30 


0.17 - 0.44 


B (a=1.80) 


0.18 - 0.23 


- 0.28 - 0.32 


0.20 - 0.46 


C (a=1.90) 


0.21 - 0.27 


- 0.30 - 0.34 


0.23 - 0.48 



^ Matching points A, B, and C with the corresponding nonlinear 
parameter a in the parenthesis (see Fig. [Jla)). 

^ Period two, two-clustered attractor (p2c2) formation occupies 
more than 90% of the basin volume in £2 — £3i and the fraction 
decreases down to 10% toward ei and £4. The rest of the basin 
volume is occupied by the formation of three or four clusters in 
random motion in ei — 62, and single-clustered chaos in £3 — £4, 
respectively. 

The prediction for the bounds based on p2c2 foliation curves. 
The upper (lower) bound from pi — f p2 (p2 — > p4) bifurcation 
threshold of the element map. 



2. Case of the Duffing GCFL 

Comparing Fig.|3]and Fig.O one clearly observes that, 
as is the case for GCML, the Duffing GCFL reduces the 
system nonlinearity caused by the averaging interaction 
via the mean field. The e ranges for the two-clustered 
regime at A,B, and C are listed in Table. HIl Again, with 
increasing nonlinearity, the e necessary to maintain the 
two-clustered attractor shifts toward the higher end. 

Let us examine how the observed parameter range in 
GCFL compares with that in GCML. 

Since two steps in GCML correspond to one cycle (27r) 
in GCFL (see (fH|) ). one may roughly predict the range 
of e in GCFL via the data of e in GCFL by a relation 

/-I cstimation\2-7r/Ai f-\ ^ \^ /'00^ 

each side expressing the contraction rate in respective 



TABLE IL The two-clustered regime in The DuffingGCFL 



Nonlinearity^ 


observed 


£z3-ran^ 




Estimation from GCMLL^ 




£1 - £2 


- £3 - 


£4 




A (fc = 0.241) 


9.3 - 10.4 


- 13.4 - 


19.0 


5.5 - 6.9 - 9.1 - 11.3 


B (fc = 0.228) 


10.1 - 11.4 


- 14.0 - 


19.5 


6.5 - 8.4 - 10.2 - 12.2 


C (fc = 0.215) 


10.8 - 12.4 


- 14.3 - 


20.1 


7.4 - 9.9 - 11.5 - 13.0 



^ The matching points A, B, and C with the corresponding 
friction parameter k in the parenthesis (Fig. [Jib)). 

^ £0 in unit of 10~^ when At = 10~^. Period one, two-clustered 
attractor (plc2) formation occupies more than 90% of the basin 
volume in £2 — £3, which decreases down to 10% toward £1 and 
£4. In £1 — £2, the rest of the basin volume is occupied by two 
clusters, each in random motion, while in £3 — £4, by a 
single-clustered chaotic attractor whose topology is chiral for 
point A and self-dual for points B and C (see insets in Fig. [Sjl. 
Estimation for ££)-range via I I22II from £inap-range in Table |l] 

models. Note that the matching of k and a is consid- 
ered between a free flow and a free map and then one 
period (27r) of a flow naturally corresponds to one itera- 
tion step of a map. On the other hand, ((22)) is considered 
between flows and maps in actual interaction, realizing 
the correspondence p^ . Tabic [Til shows that the predic- 
tion for the two-clustered attractor range roughly agrees 
with the observed ed range, although generally, the latter 
turns out to be somewhat higher than the former simple 
prediction. 

To analyze this issue globally in quantitative terms, 
we show in the right panel of Fig. [9l the band of two- 
clustered dynamics on the (A:,£d) plane, where the Ed 
axis is adjusted using the relation to facilitate the 
comparison with the GCML case in the left panel. We 
find that the two-cluster band in the Duffing GCFL runs 
in the {k,£o) plane roughly like the foliation band of p2c2 
attractor in GCML in the (a,£,nap) plane but with two 
notable differences: 

(1) the GCFL band turns out with a somewhat higher 
coupling, 

(2) the GCFL band is less tihed than the GCML one. 
These differences can be understood as follows. For sim- 
plicity, let us consider the basic two-clustered attractor 
with ■& = 0.5, which is free from iJ-bifurcation. In GCML, 
this is period two (p2c2), while in the Duffing GCFL, it 
is period one, as discussed in IIII Al (see especially (fT4)) ). 
Pairwise clusters in GCFL (each period one) emulate the 
GCML p2c2 attractor. Therefore, while the nonlinearity 
is reduced from chaos to period two in GCML. it should 
be reduced deeply down to period one in GCFL, which 
requires higher coupling at the equivalent nonlinearity. 



We are aware that this is only a crude approximation because 
all the 27r/At of the one-step contraction (1 — eu) factors are 
gathered and multiplied into a single factor after (unjustifiably) 
commuting them through all the small time (At) evolution stops. 
However, we have checked, varying At from 10^'^ down to 10^^, 
that numerically measured e o keeps the left hand side of ((22} al- 
most invariant. This supports the above approximation to some 
extent. 
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FIG. 10. The upper diagram exhibits the extended band of the two-clustered dynamics and the lower diagram exhibits full 
bifurcation trees of the elements [GCML (left) and GCFL (right)]. The format for the bands is the same as that in Fig. [9]except 
that the extended regions where the elements are periodic are depicted by hatched lines. While Fig. |9] exhibits the nonlinearity 
reduction of the chaotic elements to the periodic two-clustered attractor, the extension in this figure exhibits that of periodic 
elements to the periodic two-clustered attractor. The compositions are as follows. GCML: two-clustered attractor with period 
two (p2c2) in the central band, p2cl in Z'^, p4c4 in Z'^. GCFL: two-clustered attractor with period one (plc2) in the central 
band, p4cl,p2cl,plcl in Zg'^''''^'''^^ respectively, and p2c2 in ^4. The phases in the extended region can be easily understood 
from the bottom bifurcation diagram for the elements. For a single Duffing flow with A = 7.5, the periodic region is fc > 0.277, 
the period-two attractor {A'-_^\ A'-^) ink € (0.299,0.38), and the period-one attractor {A+, A.) in k € (0.38,0.66). 



To quantitatively verify this explanation, we have ex- 
tended the range of k (a) (Fig. [T0| such that it covers 
the periodic regime of the element flow (map). In the 
extended region, the elements are already periodic and 
hence, if the coupling is sufficiently strong to maintain 
overall coherence, a single clustered attractor is formed, 
which evolves with the same periodicity as that of the 
element {Z[ and ^^(4),(2),(i) qCML and GCFL, re- 
spectively.) If, on the other hand, the coupling becomes 
sufficiently small, a multiclustered periodic attractor is 
formed in GCML {Z2), while in GCFL, the periodic ex- 



ternal force compels the flows to synchronize in two clus- 

(2) (2) 
ters, one subject to and the other to A_ {A^) 

depending on the bifurcation state of the clement flows 
{Z^ , Z'^). In the intermediate coupling, the extension of 
the central band of Fig. IHlis formed for both GCML and 
GCFL. We clearly observe that the two-cluster band of 
GCFL is smoothly connected to the period-one window 
of the element flow at e — >■ 0, just as the GCML two- 
clustered band flows along foliation curves to the period- 
two window. 
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V. CONCLUSION 

We have constructed a Duffing GCFL as a natural ex- 
tension of the logistic GCML and investigated to what 
extent the intriguing features of the logistic GCML are 
inherited by our model. This, in turn, amounts to an 
explanation of the extended universality at the level of 
the complex systems. To compare GCFL with GCML, 
we have worked in a scheme where the element flow is 
adjusted to have the same nonlinearity as that of the el- 
ement map, in the sense that their locations in the (uni- 
versal) bifurcation trees match each other, while the cou- 
pling Ed is varied as a free parameter to seek if GCFL 
realizes the corresponding phenomena in GCML. As a 
byproduct, we have observed that the bifurcation-tree 
matching leads to remarkable agreement between the 
Lyapunov exponents of the elements over the entire non- 
linearity parameter value range. 

The phase diagram presented in Fig. [5] clearly exhibits 
that in the Duffing GCFL, the same phases as those in 
the logistic GCML are realized in the same order with 
increasing e. Furthermore, the extensive analysis of the 
Duffing two-clustered dynamics shows that the Duffing 
GCFL has the ability of iJ-bifurcation just as GCML does 
(Fig. |5]). We believe that this is the ffist systematic ob- 
servation of i?-bifurcation in coupled flow systems. This 
clearly shows that the universality between a map and 
a flow is basically extended to the level of universality 
between the systems. 

However, we have found that there is an important dis- 
tinction between the two models that comes from the fact 
that a Duffing flow admits a pair of possible attractors 



under the symmetry As a consequence, the GCFL 
attractor in the two-clustered regime consists of pairwise 
clusters (Fig. |6|) under a symmetry relation (fTTj) and each 
in period one (returns the same phase-space point at ev- 
ery 27T, modulo iJ-bifurcation). In the logistic GCML, on 
the other hand, the element map has a unique attractor 
and the two-clustered attractor consists of two clusters, 
each evolving in period two but mutually shifted by one 
step (i.e., in opposite phase). This distinction depicted 
in Fig. [7] means that the nonlinearity reduction in GCML 
is only down to period two, while in GCFL, it is reduced 
deeply down to period one. To verify this, we have nu- 
merically measured (extending the element nonlinearity 
from chaotic to periodic regime) how the band of the two- 
clustered dynamics in GCFL runs in the {k — e) space. It 
is clearly observed (Fig. [TU)) that the GCFL band flows 
into the period-one window, just as the GCML band flows 
into the period-two window. To summarize, the Duffing 
GCFL basically shares the same phases and the same 
iJ— bifurcation property with the logistic GCML; how- 
ever, detailed examination has revealed that there is a 
distinction because of the difference in the basin struc- 
tures of the elements. Finally, we are now investigating 
a GCFL of three-dimensional flows in detailF°l 

ACKNOWLEDGMENTS 

TS thanks his former student, Hayato Fujigaki, for his 
enthusiasm for this work. He and TS jointly observed the 
population bifurcation in GCFL [g^l . It has taken almost 
10 years to clarify the Duffing GCFL attractor structures 
resulting from the symmetry in the Duffing equation. 



[1] K. Kaneko, Phys. Rev. Lett. 63, 219 (1989). 

[2] K. Kaneko, Phys. Rev. Lett. 65, 1391 (1990). 

[3] K. Kaneko, Physica D41, 137 (1990). 

[4] Y. Ueda, Steady Motions Exhibited by DufSng's Equa- 
tion : A Picture Book of Regular and Chaotic Motions, 
National Institute for Fusion Science, Research report 
434, 1-12, 1980. See eq.(l) and the Duffing phase di- 
agram Fig. 1. These are also accessible in The road to 
chaos, Aerial Press (1992); see eq. (17) in page 208 and 
Fig. 7 in page 210. 

[5] M. J. Feigenbaum, Los Alamos Science 1, 4 (1980). 

[6] P. Cvitanovic in Universality in Chaos, A beautiful ac- 
count with fine figures is given in @]. 2nd edition, In- 
stitute of Physics Publishing, Bristol and Philadelphia 
(1996). 

[7] M. J. Feigenbaum, J. Stat. Phys. 21, 669 (1979). 
[8] E. N. Lorenz, Journal of the Atmospheric Sciences 20 
130 (1963). 



See a brief report of preliminary results in SecU 



[9] C. Sparrow, The Lorenz Equations: Bifurcations, Chaos, 
and Strange Attractors, New York, Springer- Verlag 
(1982). See page 99. 
[10] E. Atlee Jackson, Perspectives of nonlinear dynamics, 1, 
Cambridge University Press (1989), Fig. 7.54. See also 
page 99 in f^. 

[11] T. Moriya, T. Shimada, H. Fujigaki, Artificial Life and 

Robotics, 13, 214 (2008). 
[12] M.G. Rosenblum, A. S. Pikovsky, and J. Kurths, Phys. 

Rev. Lett. 76, 1804 (1996). 
[13] A. Pikovsky, M. Rosenblum, J. Kurths, Synchronization; 

a universal concept m nonlinear sciences, Cambridge, 

2002. 

[14] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 64, 821 
(1990). 

T. L. Carroll and L. M. Pecora, Physica 67, 126 (1993). 
[15] J. M. Kowalski, G. L. Albert, and G. W. Gross, Phys. 

Rev. A42, 6260, 1990. 
[16] H. Fujigaki, M. Nishi, and T. Shimada, Phys. Rev. E53, 

3192 (1996); H. Fujigaki and T. Shimada, Phys. Rev. 

E55, 2426 (1997). 
[17] Tokuzo Shimada and Kengo Kikuchi, Phys. Rev. E63, 

3489 (2000). 



15 



[18] T. Shimada, K. Kubo, T. Moriya, Artificial Life and 

Robotics, 14, 562 (2009). 
[19] T. Sliibata and K. Kaneko, Pfiysica D 124, 177 (1998). 
[20] T. Sliimada and S. Tsukada, Proc. of the Sixth Int. Symp. 

on Artificial Life and Robotics (AROB 6) 1 242, 2001. 
[21] W. Just, J. Stat. Phys. 79, 429, 1995. 7492 (1992). 



[22] H. Fujigaki and T. Shimada, unpubUshed, 1997; H. Fu- 
jigaki, M. Sc. thesis. School of Science and Technology, 
Meiji University, Kanagawa, Japan. 

[23] G. Perez and H. A. Cerdeira, Phys. Rev. A46, 7492 
(1992). 



